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ABSTRACT 


This program computes conventional performance parameters at required com- 
binations of speeds and flow rates. Necessary input information consists of flow areas, 
diameters, blade angles, and an estimate of design point performance. The flow 
equations and corresponding computer listings are given. A set of example calculations 
is included. 



A FORTRAN IV PROGRAM TO ESTIMATE THE OFF-DESIGN 
PERFORMANCE OF RADIAL-INFLOW TURBINES 
by Carroll A. Todd and Samuel M. Futral, Jr. 

Lewis Research Center 

SUMMARY 

A FORTRAN IV computer program for off -design performance of a radial -inflow 
turbine is presented. The thermodynamic equations used and the corresponding computer 
listings are given. Use of the program requires as input information the turbine flow 
areas, diameters, and blade angles. An estimate of design point performance is also 
necessary. The output consists of conventional performance parameters at specified 
flow conditions and speeds. A complete set of calculations for an example turbine is 
included. 


INTRODUCTION 

A procedure for predicting the off -design performance of radial-inflow turbines was 
described by Futral and Wasserbauer in reference 1. It is possible to use this proce- 
dure for hand calculations, but it is lengthy and contains several iterations. For this 
reason a digital computer program analog has been used. The computer program used 
at Lewis Research Center is given here, with the explanations necessary for its use. 

The thermodynamic equations are also given. These are the equations of reference 1, 
except for an improvement in the calculation of rotor incidence loss and the calculation 
of additional performance parameters. 

The program has been used in two ways. (1) For an existing turbine, the perform- 
ance at design operating point may be known, while information on off -design performance 
is lacking. In this case, the program provides a convenient way to estimate off -design 
performance without making actual tests. 

(2) The program is also useful as a design guide. For a proposed turbine design, 
the design point performance can usually be closely estimated as part of the design 
process. However, if the off -design operation is of sufficient importance, modifications 
in the design may be considered. In this case, the effect of a series as small design 



changes on performance may be studied by repeated use of the program. 

In addition to the engineering equations and the FORTRAN IV program, this report 
includes the input and output listings for an example turbine problem. 


ENGINEERING ANALYSIS 

The analysis consists of a one -dimensional solution of flow conditions along the mean 
streamline in the turbine, using perfect gas relations throughout. Equations were written 
in a step-by-step fashion, beginning at the turbine inlet. The analysis is written for sub- 
sonic flow only since stator choking is not expected. Symbols used in the engineering 
equations are given in appendix A, and the equations are listed in appendix B. Figure 1 
shows the location of station subscripts in the turbine. 




Figure 1. - Turbine stator and rotor. 
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Required Information 


The geometric features of the turbine which must be specified consist of flow areas, 
blade angles, and major dimensions. The flow areas needed are those at the turbine 
inlet, stator exit, rotor inlet, and rotor exit. The turbine inlet area is normal to the 
flow, upstream of the scroll. The stator exit area and the rotor entrance area are 
surfaces of circular cylinders, with elements parallel to the turbine axis. The rotor 
exit area is an annular area, perpendicular to the turbine axis. The stator exit area 
is located just inside the blade row and the rotor inlet and exit areas are calculated at 
locations just outside the blade rows. 

The blade angles at the stator exit and rotor exit must also be specified. These are 
angles at the mean streamline. 

Diameters needed are the stator exit diameter, rotor inlet tip diameter, and the 
mean diameter at the rotor exit. The axial blade length at the rotor tip is also neces- 
sary. 

The required inlet conditions of the working fluid are the total temperature and total 
pressure at the turbine inlet. The gas properties needed are the specific gas constant, 
the specific heat at constant pressure, and the specific heat ratio. 

Input information which applies to the example turbine is as follows: 


Turbine inlet area, sq ft 0. 12989 

Stator inlet area, sq ft 0. 04078 

Rotor entrance area, sq ft 0. 04175 

Rotor exit area, sq ft 0. 04719 

Stator exit diameter, sq ft 0. 42650 

Rotor inlet diameter, sq ft 0. 41417 

Stator exit axial dimension, ft 0. 03208 

Stator exit angle, deg 72.47 

Rotor exit angle, deg -56. 85 

Number of blades at rotor tip 22 

Turbine inlet total pressure, °R 2060 

The thermodynamic properties of the working fluid are as follows: 

Specific gas constant 18. 437 

Specific heat ratio 5/3 

Constant pressure specific heat 0. 05924 
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The design point quantities, used in determination of frictional coefficients, are the 
following: 


Mass flow rate, lb/ sec 0. 7484 

Total efficiency 0. 897 

Wheel speed, rpm 36 000 

Total pressure ratio across stator 0. 98 

Total pressure ratio across turbine 1. 7397 


Operation of the Program 

In solving the equations in appendix B, the computer uses the equations in essenti- 
ally the order listed. Solutions for the quantities ^V/V^q, and ^V/V c ^ ^ 

are accomplished by means of the Newton-Ralphson method since these quantites are 
implicitly defined. 

For operation of the program, the wheel speeds of interest are specified in advance. 
The initial value for (v/V cr ) p the velocity ratio at the stator exit, as well as its incre- 
ment size are also specified. The final value for this quantity will not be known and is 
determined by the computer for each wheel speed. For any given wheel speed, repeated 
increases in the value of fv/V cr j ^ will eventually result in a flow rate which the turbine 
cannot pass. Mathematically, this condition is evident in that no solution exists for 
^W/W cr j^, the relative velocity ratio at the rotor exit. When this occurs, the maximum 
value for y/V c ^J ^ is computed to five significant figures. The computer prints out a 
complete set of answers for each value of (v/V cr ^ j for any given wheel speed. When 
one wheel speed is finished, the computer goes on to the next speed until all specified 
speeds have been used. 

In this discussion, it is always assumed that the temperature and pressure of the 
working fluid entering the turbine are constant. This has fulfilled all needs thus far. 
However, if the program should be required to operate in some other manner, the 
amount of additional instructions would not be great. 


Loss Coefficients 

Calculation of rotor incidence loss makes use of the angle cp in equation (B17). 

This is the optimum incidence angle and is that incidence angle associated with minimum 
rotor incidence loss. The calculation of cp follows the method of Stanitz (ref. 2), 
developed to account for ’’slip" in centrifugal compressors having radial blading. 
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The loss coefficient kj applies to the stator and must be evaluated before the main 
program is used. The determination of kj^ involves simultaneous solution of equations 
(Bl) to (B5) using the design value of mass flow rate. For this purpose, the total pres- 
sure ratio across the stator Pj/Pq is determined from design data. A separate routine 
is provided for this solution. 

The determination of the loss coefficients m and k were presented as a graphical 
procedure in reference 1. In the computer determination, this is done by means of a 
search routine. By repeated calculations, the computer finds a pair of values for m 
and k which yield design values for mass flow rate, efficiency, and pressure ratio. 
Design speed and inlet conditions are specified. 


Output Information 

A sample of the computer output is shown for the example turbine. To assist in 
reading this, a list of machine output symbols in terms of engineering symbols is 
included in the section on program output. The information contained in the machine 
output may be classified under several headings as follows: 

(1) Input information - For each speed, the following quantities are printed out: 

(a) Design speed 

(b) Percent design speed 

(c) Fluid inlet temperature and pressure 

(d) Turbine dimensions and angles 

(e) Fluid constants 

(f) Number of rotor blades (at tip) 

(2) Performance parameters 

(a) Efficiences 

(b) Blade -jet speed ratio 

(c) Pressure ratios 

(d) Specific work 

(e) Ideal work 

(3) Internal variables - Internal variables are listed for each station throughout the 

turbine as follows: 

(a) Pressures 

(b) Temperatures 

(c) Velocity components 

(4) Additional quantities - Certain ratios have been found useful and are included 

in the print -out, such as the following: 
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(a) wNe/6 

(b) N/v^ 

(c) w/TJ/p^ 

(d) r/p' n 

(e) Losses calculated on a pound mass basis denoted by L with appropriate 

subscript 

(f) Fractional losses denoted by A 77 where the loss has been divided by the 

ideal work 

(g) Losses in efficiency, given in equation (B 31 ) as r^, rj R , and rj j, which 

are useful in studying loss distributions in the turbine 


THE FORTRAN IV PROGRAM 

In this section are presented the FORTRAN IV routines, input and output consid- 
erations, and the internal documentation. The program operates in three distinct 
modes, each depending on the type of information available or the desired output. In 
mode one operation, the loss coefficients are assumed known. Off-design calculations 
are made by means of changes in speed and changes in the value of (V/V cr ^. This 
mode furnishes performance results, once the loss coefficients are known. In mode 
two operation, sets of loss coefficients are entered and calculations are made at design 
values of speed and (v/V cr )i* This mode is used for studying the results of variations 
in frictional effects. In mode three, the program determines the loss coefficients and 
produces one complete set of calculations representing design conditions of the turbine. 
This set of calculations verifies the determination of the loss coefficients. 


Program Input 

Actual input for the example with card format is presented in this section, followed 
by the definitions of the terms used on the cards. 
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title 


PROJECT NUMBEI 


ANALYST 


SHEET, OF 


STATEMENT 

NUMBER 

z 

© 

FORTRAN STATEMENT 

IDENTIFICATION 

1 2 5 4 5 

6 

7 8 9 10 It 12 

— 
13 14 IS 16 17 18 

19 20 21 22 23 24 

ZS 26 27 28 29 sd 31 32 33 34 3S 36 

37 38 39 40 AI U2 43 44 4S 46 47 4B 49 SO 51 52 S3 54 55 56 57 58 59 60 

61 62 63 64 65 66 67 68 69 70 71 72 

73 74 75 76 77 78 79 81 






MODE ONE 

FORMAT 





















PREF 


WREF 











.98 66 

7 

. 748 











AO 


A 1 


' A V 


A 4 

D 1 


D.3. 


XL . 

. AL 

1, , 

.05 54 

8 

. 04078 

.0 417,5 

.047188 .4265 . 

.4.141 6 6 . 

. 0.3.2.08 3. . . .7.2 


B 4 


' ZBR 

PO 


T 0 

U 3 


U4 


G 

] R 


-56.5 

5 

22 


358 8 

. 12 

2 060. 

780.6838 

443 . 

245 

1 . 6667 

, , 18 

.43 7, ( 

CP 


VST ART 

DEL V 


VEND 








.0592 

4 

2 .£ 


• °, 5 


.,9, , , , 








MODE 














» -t -f i — 
1 r . 














XM^ 4 ^ 


XK 










2.9 21 

4 

01 .0950602 











PER.C.T 














120^1 




no. , 




10,0. , 




_ 8 ( 0^_ 





/V 
















ion 








ip^i 











' 





1 I 1 t 1 




















m 

i 

7 0 9 10 II 12 

13 14 IS 16 17 18 

19 20 21 22 23 24 

25 26 27 28 29 30 

31 32 33 34 35 36 

37 38 39 40 41 42 

43 44 45 46 47 48 

49 50 51 52 53 54 

55 56 57 58 59 60 

61 62 63 64 65 66 

67 68 69 70 71 72 

73 74 75 76 77 78 79 • 


NASA-C-836 (REV. 9-14-59) 




















The definitions of the terms used on the cards are as follows: 


Card 1: 
PREF 
WREF 
Card 2: 

AO 
A1 
A3 
A4 
D1 
D3 
XL 
AL1 
Card 3: 

B4 

ZBR 

PO 

TO 

U3 

U4 

G 

R 

Card 4: 

CP 

VSTART 

DELV 

VEND 


estimate of Pj/Pq to determine 
estimate of W to determine k^ 

area of turbine inlet, A q , sq ft 
area upstream of stator exit, Aj, sq ft 
area upstream of rotor inlet, Ag, sq ft 
area downstream of stator exit, A 4 , sq ft 
diameter upstream of stator exit, D^, ft 
diameter upstream of rotor inlet, Dg, ft 
axial distance, l , ft 
absolute gas angle, a^, deg 

free stream gas flow angle, /3^, deg 
number of rotor blades at tip of rotor, Z 
inlet pressure, p^, lb/ sq ft 
inlet temperature, T^, °R 

design speed upstream of rotor inlet, Ug, ft/ sec 
design speed downstream of rotor exit, U^, ft/ sec 

ratio of specific heat at constant pressure to specific heat at 
constant volume, y 

specific gas constant, R, ft-lb/ (lb) (°R) 

specific heat at constant pressure, c , Btu/(lb)(°R) 
starting value of (v/V cr ^ 
incremental value of (y/V cr )i 
end value of (v/V cr )i 
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1 1 


Card 5: 
MODE 


Card 6: 

XM 
XK 
Card 7+: 
PERCT 
Card 6+: 

XM 
XK 
Card 6: 

PDSGN 

ETAD 

P 

E 


if MODE = 1, 
if MODE = 2, 
if MODE = 3, 


a performance run varying percent speed 
manual determination of m and k 
automatic determination of m and k 


loss coefficient 

loss coefficients > 


percent of design speed J 

-\ 

loss coefficient > 

loss coefficient 

J 


specified P0/P4 or pjj/p^ 
specified or ri s 
if p = 1, Pq/p 4 is specified 
if p = 2, Pq/p 4 is specified 
if E = 1, 77 D is specified 

'S 

if E = 2, Tjj. is specified 




MODE ONE 


MODE TWO 


MODE THREE 


Description of the Output 


A list of the machine output symbols, internal program variables, and their corre- 
sponding engineering symbols is given here. Refer to appendix A for the physical 
interpretation. Following this are sample output listings showing the results of 
calculations for modes one, two, and three using the input data given previously. It 
should be noted that each block of output is for a specific value of (v/V cr ^. 
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Machine Output Symbols and Definitions 


ALPHA2 

“2 

BETA3 

% 

DEE 

a? ?e 

DEI 

Arjj 

DELTA 

6 

DER 

At?r 

DES 

A77 s 

EQ -SPEED 

N eq 

ETA-E 

VE 

ETA-I 

n 

ETAR 

% 

ETA-S 

% 

ETA-T 


LE 

l e 

LI 

L I 

NOT 

n/VTS 

NU 

V 

PO,/P4 

P0/P4 

PO,/P4, 

Pb/P ’4 
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Pl,/PO, 

P2,/PO, 

P2,/P1, 

P4/PO, 

P4/P4ID 

P4/P4, 

P4./PO, 

PHI 

PV/PVCR)0 

PV/PVCR)1 

PV/PVCR)2 

PVX/PVCR)3 

PW/PWCR)4 

R2 

RH0)0 

RH0-VCR)0 

ROTOR-LOSS 

SMALL -W 
SPEED 
SPEED -PAR 


p’i/po 

P^P’o 

P2/P1 

P4/P0 

PyP’4,id 

P4/P4 

P 4 / P’o 

<P 

(pv/p’v cr ) 0 
(pv/p'v cr ) ! 
(pV/p’V cr ) 2 
(pV x /p’V C r ) 3 
(pw/p"w cr ) 4 

r 2 

f>0 

(P' v cr )0 

w 

N 

wN/ 6 



STATOR-LOSS 

L s 

T, ,/T, )3 


T3, , 

rptt 

1 3 

T4, 

T 4 

T4,/T4,, 

(t’/t ") 4 

T4,,/T3,, 

T 4 /T 3 

THETA-CR 

®cr 

TOP 

r/Po 

V3 

V 3 

V4 

V 4 

V/VCR)0 

( V / V cr) 0 

V/VCR)1 

( V/V cr)l 

V/VCR)2 

( V / V cr) 2 

V/VCR)3 

( V / V cr) 3 

VCR)0 

V cr,0 

VCR4 

V cr,4 

VU3 

CO 

s/ 

VU4 

V u,4 

VU/VCR)3 

(V V cr) 

VX3 

V x,3 
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VX/VCR)3 

(V V cr) 3 

W1 

Ah’ 

W3 

W 3 

W4 

W 4 

W/WCR)3 

( W / W cr) 3 

W/WCR)4 

( W /W cr ) 4 

WCR3 

W cr,3 

WCR4 

W cr,4 

WEQ 

w eq 

WID 

Ah id, s 

W, ID 

Ah id,t 

WTOP 

wVTq/Pq 



Mode One Sample Output 


OFF DESIGN PERFORMANCE OF A RADIAL INFLOW TURBINE 


TURBINE GEOMETRY 

AO 

0. 12Sge7 

A1 

0. 407800 E—0 1 

A3 

0. 417450 E-01 

A4 

0.47 1880E— 0 l 

Dl 

0.42 6500 


D3 

0.414166 

L 

0. 32 0830 E— 01 

AL1 

72.47000 

B4 

-56.85000 

Z 

22.00000 

INLET CCNDIT ICNS 

PO 

36 CO. 000 

TO 

2060.000 







DESIGN SPEEDS 

U3 

780.6838 

U4 

443.2450 







FLUID PROPERTIES 

G 

1.666700 

R 

18.43700 

CP 

0.592420 E-01 

EPS 

0.910949 




MODE 1 


CES IGN 
CESI GN 


PERFORMANCE RIJN VARYING PERCENT SPEEO 
M 3.592625 
K Q.754039E-01 
Kl 0. 763553E— 0 I 

PERCENT lOO.COOC 
IJ 3 780.6638 

U4 443.2450 
PW/PO, 0.9BC000 
SMALL-W 0.748400 


V/VCR ) 1 

0.2COOOO 

P2 */PQ* 

0.998079 

V 3 

253.5031 

WCR 3 

1258.994 

VCR JO 

1235.911 

P2./P1 , 

1 .000023 

BETA3 

-82.29354 

W4 

132. 7295 

Rhfl jn 

0 • 94 7862E— 0 1 

ALPHA2 

73.31849 

T » t/T» ) 3 

L. 037702 

WCR 4 

1217.305 

RHO-VCR )0 

1 17. 1 A 7 3 

R2 

0.212270 

T3, , 

2137.666 

VU4 

332.1319 

PV /P VCR ) 1 

0.1S70C8 

VU/VCR ) 3 

0.196464 

T4../T3, , 

0.934870 

T4./T4, , 

1.016528 

SMALL-W 

0. 283144 

VX/ VCR J 3 

0.5 89386E-01 

W/WCR) 3 

0.431123 

P4 f /P0, 0 

.885221 

PV /P VCR JO 

0.1 86084E-01 

V/VCRJ3 

0.205114 

PVX/PVCR) 3 

0.580112 E-01 

VCR4 

1227.323 

V/VfR JO 

C . 1 (6 108E-01 

VU3 

242.8122 

P4/P4I D 

0. 918390 

V4 

339. 9744 

V/VCR ) 2 

C.20C089 

W3 

542.7817 

PW/PWCR) 4 

0.108550 

P4/P4* 

0. 952730 

PV/PVCRJ2 

0.197093 

VX3 

72. 84285 

W/WCRJ4 

0.109036 

P4/P0* 

0. R43376 

W1 

1.651284 

ETA-E 

0.497547 

W, ID 

5.808875 

T4» 

2031.472 

WID 

E.C38581 

ETA-T 

0.291155 

ETA-S 

0.210396 

WEQ 

0.183813 

ETA- I 

0.778959 

LE 

2.308287 

DEE 

0.287151 

LI 

2.267151 

ROTOR-LOSS 

1.68918? 

STATOR-LOSS 

0 .92886 8E—0 1 

DER 

0.210134 

DES 

0. 1 1 555 1 E-01 

SPEED 

36000.00 

EO-SPEED 

29695.14 

MU 

1.730507 

THET4-CR 

1.469718 

DELTA 

1.7CI 146 

SPEED-PAR 

5991.953 

Wl-EQ 

1.150754 

P0./P4, 

1.12966? 

NOT 

793.1747 

WT OP 

0.514045 

TOP 

0 . 474460E- 01 

DEI 

0.281412 

P0./P4 

1. 18571 0 

PHI 

-43.96986 

PU/PO, 

0.998056 

ETA-R 

1.0661 10 

NSP 

12C .8699 







V/VfR) l 

0 . 25000 C 

P2 ./PO. 

0.996982 

V 3 

316.8903 

WCR 3 

1249.548 

VCR JO 

1235.9 1 1 

P2./P1. 

1.000035 

BET A 3 

-79.19788 

W4 

167.5706 

RHOJO 

0.947862E-01 

ALPHA2 

73.31857 

T, ,/T, ) 3 

1.027189 

WCR4 

1207.532 

RHO-VCR JO 

1 1 7. 1473 

R2 

0.212770 

T 3f , 

2105.710 

VU4 

302.9650 

PV/PVCR J I 

C. 244164 

VU/VCR) 3 

0.245580 

T4» * /T3» , 

0.933882 

T4* /T4f » 

1.012364 

SMALL-W 

0.350528 

VX /VCR ) 3 

0.7 37055E— 0 1 

W/WCR) 3 

0.388769 

P4 f /PD, 0. 

847779 

PV/PVfR JO 

C . 230369E- 0 1 

V/VCR) 3 

0.256402 

PVX/PVCR) 3 

0.7189 59 E-01 

VCR4 

1214. 974 

V/VCR JO 

0.2 304 1 5E- 0 l 

VU3 

303.5152 

P4/P4I D 

0.926179 

V4 

316.5768 

V/VCR J2 

C. 25011 1 

w3 

485.7858 

PW/PWCR ) 4 

0.137770 

P4/P4, 

0.95811 8 

PV/PVCR )? 

0.244267 

VX3 

91.09343 

W/WCRJ4 

0.138771 

P4/P0, 

0.812273 

W 1 

4. 1C0490 

ETA-E 

0.626467 

W* I D 

7.800876 

T4, 

1990.797 

WIO 

5.739316 

ETA-T 

0.525645 

ETA-S 

0.421024 

WFQ 

0.277558 

ETA-I 

C. 839335 

LE 

7.000868 

DEE 

0.205442 

LI 

2.073232 

ROTOR-LOSS 

1 .428626 

STATOR-LOSS 

0. 145112 

DER 

0.146686 

DES 

0.148996E-01 

SPTED 

36CC0.00 

60-SPEED 

29695. 14 

NU 

1. 117917 

THETA-CR 

1.469718 

DELTA 

1.7CU46 

SPEED-PAR 

7417.945 

Wl-EO 

2.789984 

PD./P4, 

1.179552 

NOT 

793.1747 

WTOP 

0.636380 

TOP 

0 . 142408 

DEI 

0.212872 

P0./P4 

1.231113 

PHI 

-37.64624 

PI ./PO, 

0.996947 

ETA-R 

1.044782 
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Mode Two Sample Output 


OFF DESIGN PERFURMANCE QF A RADIAL INFLOW TURBINE 


TURBINE uEUMETRY 

INLET LuNDlT IUnS 
DESIGN SPEEDS 
FLUID PKUPeRTIcS 


AU 

0. 12 98d 7 

03 

0. 414160 

Pu 

3600. ODD 

U3 

780.6838 

0 

l.O6o700 


ai q.40730oe-oi 

L 0. 320330E-01 

TU 2060.000 

U4 443.2450 

R 18.43700 


A3 0.417450E-01 

AL1 72.47000 

CP 0 .592420 £-0 1 


A4 0.471880E- 

84 -56.85000 

EPS 0.910949 


01 0.425500 

Z 22.00000 


MODE 2* MANUAL 
M 
K 
K1 

PERCENT 

UJ 

U4 

DESIGN Pi./PU, 
DESIGN SMALL-W 


DETERMINATION UF 
3.0 UuUUO 
D .5000 DUE- 0 1 
U.76385JE-01 
lOU.UJOO 
780.ob38 
443.2450 
0. 980000 
0.748400 


M AND K 


v/vcri a 

0 • 20 OU 00 

P2 ./PU, 

0.998695 

vi 

256.3198 

WCR3 

1258.537 

VCR )U 

1233.911 

P2, /PI , 

1.000640 

BETA3 

-82.25358 

U4 

130.5060 

KHUIU 

0. 9478O2L-01 

ALPHA2 

73.51138 

T. ,/T, >3 

1.036950 

WCR4 

1216.632 

RHO-VCk 10 

11?. 14? 3 

R 2 

0.2L2D38 

T3, . 

2136.116 

VU4 

333.9932 

PV /P VCR i 1 

0 . 1 9 70 Uo 

VU/VCK) 3 

0.193847 

T4 » 1 / T 3» , 

0.934823 

T4./T4, , 

1 . 016820 

SMALL— M 

0.263144 

VX/VCK1 3 

0.589232E-0 1 

fit d CR)3 

0.428959 

3 4 » /P3 » 0. 

900118 

PV/PVLK »J 

O. 1 8 8064b — 01 

V/ VCR i 3 

0.207393 

PVX/PVCR) 3 

0 • 579754E— 01 

VCR4 

1227.023 

V/VLK )U 

0 . 1861 06 c - 01 

Vu3 

245.7571 

P4/P4I D 

0.934413 

V4 

341.5363 

V/VCK J<£ 

J. 202 D 16 

m3 

539.8610 

PW/PWCR) 4 

0.106788 

P4/P4, 

0.952Z7B 

PV/PV CK )2 

0.199409 

VX3 

72. 82386 

W/rfCRl 4 

0.107251 

P4/P0, 

0.857163 

Ml 

1. 730136 

b TA-E 

0.339111 

W » I D 

5.030349 

T4» 

2030.478 

MAD 

/ • 29o7 7o 

ETA-T 

0.34/920 

ETA-S 

0.239853 

we a 

0. 183813 

ETa-1 

0 . 863162 

LE 

2.329552 

DEE 

0.319258 

LI 

2.233187 

RoTUK-LUSo 

0 • 924096 

S f A TuR- LOSS 

0.630S30E-D1 

OER 

0.126645 

DES 

0.864533E-02 

SPeEO 

36000.00 

EU— SPEtD 

29695.14 

NJ 

1.291541 

THE TA-CR 

1.459718 

DELTA 

1. 70 1 146 

8 PEED— PAR 

3991. 933 

Wl-EO 

1.190812 

P0./P4, 

1.110965 

NOT 

793. 174/ 

WT OP 

0.514045 

TOP 

0.490976E-01 

OEI 

0.306051 

PU./P4 

NSP 

1 • lo 60 40 
133. 3119 

PHI 

-43.97732 

P1./P3, 

0.998056 

ETA-R 

1. 184420 


V/VLK ) 1 

0 • *i5 UUU U 

P2 ,/PO, 

0.997949 

V 3 

320.4066 

WCR3 

1248. 973 

VCR Id 

\.£.yj .9ii 

P2./P1, 

1.001005 

BET A3 

-79.12022 

W4 

165.0909 

K Hu ID 

0.947bO2E-01 

AlPHAE 

73.61363 

T../T, 13 

1.021249 

WCR4 

1206. 937 

KHU-VCKIJ 

117.1473 

K2 

0.212065 

T 3* t 

2103*772 

VU4 

305. 0409 

PV/PVLK 1 1 

D .*14416* 

VU/ VL*I 3 

0. 243 ? 58 

T4, f /r3, . 

0.933821 

T4./T4,, 

1.012692 

SMALL -M 

0 . J6U52 b 

VA/VCRJ 3 

0./36733E-01 

H/WCRI3 

0.386049 

P4./P3, 0. 

850185 

PV/FVLR Id 

0 • 23D869l- 01 

V/ VCR) 3 

0.259247 

PVX/PVCR) 3 

0.71 826 3E-0 1 

v;R4 

1214.572 

V/VLK Id 

0.2J0416E-D1 

VD3 

307.1957 

P4/P4ID 

0.940376 

V4 

318.1267 

V/VLK it. 

0 .253144 

*3 

482 . 1840 

Prf/PWCR) 4 

0.136826 

P4/P4, 

0. 957672 

PV/PVCKJ2 

D.24 7Do6 

VX3 

91. 03610 

W/ WCRJ 4 

0.136785 

P 4V P □ « 

0.823775 

Ml 

4.17 8 303 

E TA-E 

0.O8U8J7 

W * I D 

7.135124 

T4, 

1989.480 

n ID 

9. 103924 

E TA-T 

U. 583624 

ETA-S 

0.458877 

WE 3 

0.227558 

ETa-a 

0. *0434o 

LE 

2.021 143 

DEE 

0.221959 

LI 

2.035274 

RUT UR-l JSS 

J.7/8U77 

S/ATUK-tUSS 

D. 983 3 12 1— 0 1 

DER 

0 • 854474E-0 l 

OES 

0. 108228E-01 

SPEED 

3uDoO. DO 

EU-SPEEU 

29693.14 

NJ 

1.156144 

THET4-SR 

1.469718 

0«_L I A 

1. 701146 

SPEhU-PAR 

7417.943 

Wl-EQ 

2.843064 

P0./P4, 

1. 162541 

NDl 

793. 1747 

WTUP 

0.636380 

TOP 

0.1451L7 

OEI 

0.223511 

PU, / P 4 

1.^13924 

Pn L 

-37.65760 

Pl./PO, 

0.996947 

E T4- R 

1.126307 
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Mode Three Sample Output 


OFF DESIGN PERFORMANCE OF A RADIAL INFLOW TJRBIME 


TURBINE GFCHETRY 

AO 

0.129887 

A1 

0 • 407800E— 01 

A3 

0. 417450 E-01 

A4 

0.47 1B80E— 01 

ni 

0.426500 


D 3 

0.414166 

L 

0. 3208 30 E— 01 

AL1 

72.47000 

B4 

-56.85000 

z 

22.00000 

INLET CON HIT IONS 

PO 

3600.000 

TO 

2060.000 







DESIGN SPEEDS 

U3 

780.6838 

U4 

443.2450 







FLUID PROPERTIES 

G 

1.666700 

R 

18.43700 

CP 

0.592420 E— 0 1 

EPS 

0.910949 




MODE 3. AUTOMATIC DETERMINATION OF M AND K 
M 3. 94783 5 
K 0. 68367IE— 01 
K 1 0. 7 63553E- 01 


PERCENT 100.0000 

U3 780.6838 

IJ 4 443.2450 

DESIGN Pl./PO, 0.980000 

DESIGN SMALL— W 0.748400 

OESIGN P0./P4. 1.740000 

DESIGN FTA-T 0.913000 

V/VCRH C. 615818 

P2 t/PO, 

0.981934 

V 3 

783.3564 

WCR 3 

1177.725 

VCR in 

1235.911 

P2./P1 , 

1.001973 

BETA3 

-7.71 1356 

W4 

444.0692 

PHO 10 

0.94 7862E— 01 

ALPHA2 

73.38271 

T t * /T, ) 3 

0.908058 

WCR4 

1133.050 

RHO-VCR in 

127.1473 

R2 

0.212199 

T3t , 

1870.599 

VU4 

71.49717 

PV /P VCR 11 

0.530352 

VU/VCR) 3 

0.606932 

T4../T3, , 

0.925572 

T4,/T4, t 

0. 974083 

SMALL -W 

0.748445 

VX/VCR13 

0.1 82683 

W/WCR) 3 

0.193457 

p 4,/po, o 

.575019 

PV/PVCR in 

0 • 491 383E- 01 

V/VCR >3 

0.633829 

PVX/PVCR) 3 

0.155864 

VC 8 4 

1118.271 

v/vcr in 

0.49233 IE- 01 

VU3 

750.1139 

P4/P4I D 

0.965595 

V4 

253.2051 

V/VCR 12 

0.618131 

W3 

227.8395 

PW/PWCR) 4 

0.369567 

P4/P4, 

0.968264 

PV/PVCR 12 

0.531715 

VX3 

225.7794 

W/WCR) 4 

0.391924 

»4/PH, 

0.556771 

W1 

22. 12419 

ETA— E 

0 .9 1 83 37 

W.ID 

24.23218 

T4 * 

1686.602 

win 

25.48582 

ETA-T 

0.913009 

ETA-S 

0.868098 

W6Q 

0.485880 

FTA-I 

0.919463 

LE 

1.280390 

DEE 

0 .502393E-01 

LI 

0. 2 868 48E-01 

RHTnR— LOSS 

1. 342740 

STATOR-LOSS 

0.801912 

06 R 

0. 526858 E-01 

OES 

0.3146 50E-01 

SPEFD 

36000.00 

60-SPEED 

29695. 14 

NU 

0.691074 

THETA-CR 

1.469718 

DEL TA 

1.701146 

SPEED-PAR 

15838.75 

Wi-EO 

15.05336 

Pn,/P4, 

1*739072 ' 

NOT 

793.1747 

WTCtP 

1.358794 

TOP 

1 .640600 

DEI 

0. 1 1 2552E- 02 

PO. /P 4 

1. 796072 

PHI 

-17.2871 8 

Pl,/PO, 

0.980000 

FTA-R 

0. 969702 

NSP 

75.74157 
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Main Program Listing 


C A FORTRAN IV COMPUTER ROUTINE TO PREDICT OFF DESIGN PERFORMANCE 
C OF A RADIAL INFLOW TURBINE 
C 

SIJBCX. G)=X*( l.-(G-U)/ (G+l. )*X*XI**( l./(G-l« ) ) 2 

DATA AC/32. 1739/, P I /3. 14159/, XJ/778. 16/ ,Fi/. 01745/ 

CAT A TOL/ *000001/ , XNK/ 5. / , XMN/9. / 

COMMON /GFTK/PREF.WRFF ,G . AC , A 1 • AL l , TO , R , PO 
It POSGNtCPt GlfG2#G3*ETAD*U3»L4tVU3*B4tT4T3»WCR3 *TT3 ,W3WC3,B3, 

2 PHItPPPPIPtP IP POP 
INTEGER PTEST , ETEST 
RFADI 5 , 4 G C ) PR E F * WR F F 

INPUT TURBINE GEOMETRY- I NL E T C CND l TI ON S-SPE E DS-FLU I D PROPERTIES 

4 

1 R E ADI 5*400) AO , A 1* A 3 , A A ,D 1 .03 , XL ,A LI , B4 , ZBR *PO,TO, 

1U 3 * tJ4. G.R.CP. VSTART ,0ELV1 ,VEND 5 

li 2 1 = U 3 
U41=U4 

CALL GETKl(XKXtVlVClO) 

6 

G I = G + 1 . 

G 2 =G— 1 . 

G ‘*=G2/ G 1 

EPS=1.4/G*t (G1/2.)**IG/G2 1/1.8929) 8 

R 1 = D 1/ 2 . 

R 3- C3/ 2 . 

DETERMINE MODE OF OPERATION 

READ(5*40I) MODE 
C MG DF= 1 * IMPLIES A PERFORMANCE RUN VARYING PERCENT SPEED 
C M0DE = 2* IMPL IES INPUTTING VARIOUS VALUES OF M AND K AT DESIGN SPEED 
C MO CF= 3* IMPL IES LETTING THF MACHINE DETERMINE THE VALUE OF M AND K 


C 9 

30 Gr TO ( 2, 3.4) , MODE 

2 R F ADI 5*400) XM • XK 12 

5 READ! 5*400) PFRCT 13 

V 1 VC 1 = V S TAR T 
DEL V =DEL V 1 
GO TO 6 

3 READ! 5*400 XM *XK 15 

PERCT- 100. 

V 1VC 1 = V START 
D EL V =DE L V 1 
GO TO 6 

4 READ! 5*405) PD SGN, E TAD , PTEST ,E TE ST 17 

405 FORMAT ( 2F10.4.2I l) 

X K= I .E— 8 
K = 1 

VIVC1=V1VC10 
P ERCT = 100 . 
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n o r> on 


6 U 3 =U 31*PERCT/100# 

U4=U41#P tRCT/ 100. 

C C ASF CONSTANTS 
GO TO 40 

50 VCR0=S0KT( 2 . /G.l*AC* TG*G*R J 
RHOOP=PO/R/TO 
PVCRO=VC RO*RHOOP 
C STATOR ANALYSIS 

CriSAl l=COS(ALl*Fl) 

S INAL1 = SIN(AL1*F1) 

C 

2 3 CrNTIMUE 

RhnVCl=SUB(VlVCl.G) 

0 1= 2 . *XK 1*RH0VC 1*A1*C0SALI /A0/G1 

V = V IVC 1 
X = Y /? • 

0 1=01/ ( 1 .-Y*Y) 

10 P1PP0P=< 1 .-XK J*G3*(X*X+Y*Y)/(1.-G3*Y*Y) )**(G/G2) 

SU = P 1 PPOP*RHOVC 1*P VCRO*A 1 #COSAL 1 
R HO VCO = SW/PVC R O/AO 
F = SUB( X, C-J-RHOVCO 
0 = ( F+R KJVCO )/X 

FP = 0-2.*X**2/G1*Q**(2.-G> +X*01*P1PP0P**I2.-G) 

X I=X-F/FP 

IFURSf { X 1-X)/XU.LT.T0L) GO TO 9 
X=X I 

GO TO 10 

STATOR TO ROTOR ANALYSIS 
9 VGVC0=X1 

V 2VC 2= ( XK 1 + l« ) /( XK+1* ) *Y*Y-<XK-XK1 ) /( XK+1. )*X*X 
V2VC?=SQRT I V2VC 2 ) 

R F 2 VC 2= SIJB ( V2VC 2 *G ) 

P2PP0P=( 1 G3*XK* ( V0VC0**2+V2VC2**2) / ( 1 • -G3 *V2 VC2 **2 ) ) )** (G/G2) 
P2PP IP = P 2PP0P/P 1PP0P 

TANAL?=2.*P I*R l*XL*RH2 VC2*P2PP 1P*SINAL1/ (RHCVC l*Al*CGS AL1) 

A L 2 = A TAN ( T ANAL 2 ) 

SINAL2=SIN( AL2) 

R 2=R l* S I NAL 1/S INAL2 

ROTOR INI. ET ANALYSIS 

VLVC 3=R 2/R 3*V2VC2*S I NAL2 
X=VlVCl/2. 

12 0=R?/;<3*RH2VC2*C0S1AL2> 

F=0*( l.-G3*<X*X+VUVC 3 ** 2 ) )**C-1./G 2)-X 
FP*2./G1 *X*Q**2/CF+X)- 1. 

X 1=X-F/FP 

IFCABSC <Xi-X)/Xl).LT.TCL> GO TO 11 
X = X1 

GC TO 12 
11 VX VC3= X 1 

V3VC3-SORT (VXVC 3**2+ VUVC 3**2) 

VU 3=VC RO *VUVC3 
VX3 = VCRfl*VXVC3 

V 3= VCRO* V3VC 3 


numamikvii ■■ ■■■■!■■ ^^MIBIIIIII HIM 


I HI HIM I ■! liniUlllllHIH 
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C 


C 

C 


r 

c 

c 


W 2 = < V3**2-2„*U3*VU3+U3**2 )**. 5 
CDSB3= VX 3/ W 3 
S INB3= ( VU3— U3) /W3 

B 3 = AR S IN { S I NB 3 ) 51 

T 73 = 1 .-G3*{ 2.*U3*VU3-U3**2 ) /VC R0**2 
T 3PP=TO*TT3 

T4T3=l.-( U3**2-U4**2) / ( 2 .* AC* X J*C P*T3 P P ) 

W3WC3=VXVC3/( SQRT( TT3 ) *C0 SB3 ) 52 

WCR3=W3/W3WC3 

RFXVC3-SUB( V3VC3*G ) * VX VC 3/V3VC3 

WX3=W3*COSB3 

V U30PT =U 3* ( 1.- 1.98/ZBR) 

W U 30P7 = V li 30PT— U3 
PHI=ATAN(fcU30P7/WX3) 

53 

I F l MODE ,NF *3 )G0 TO 300 

CALL GFTMi XK*XM tltETAD.l) 57 

300 CONTINUE 
ROTOR AN Al YSIS 


X=V1VC 1/2. 

0=RHXVC3*A3/( A4*C0S(B4*F 1) ) / ( T T3 * T4 T3 ) ** ( G1 /G2 / 2 . ) 59 

0 1 = W3*IC3**2*< XM*XK+SIN(6 3- PHI )**2) /T4T3 61 

K CWN T= 0 

15 Q2=1.-G3*X*X 
KCV*NT=KOVsNT +1 
Q3=XM*XK*X*X+Q 1 

F 2 = C l G3* ( Q 3 )/Q2)**<G/G2> 64 

F 3 = SUB ( X * G ) 

F =0 / F2— F 3 

F 2P =G/G2*F 2** ( 1 ./G)* ( - 2 .*G 3**2 *03/ Q2 **2*X-2. *G3 *XM*XK* X/ Q2 ) 65 

F3P=-2.*X**7/G1*(F3/X)**( 2.-GI+F3/ X 66 

FP=-Q/F2**2*F2P-F3P 
X 1 = X— F/FP 


CHFCK FOR CHOKING CONDITIONS 

IF( ( XL .GT . 1. ) .OR. ( XI. L T.O. ) .OR. ( KOWNT • GT. 300 I ) GO TO 13 

IF(ABSl (Xl-X)/Xl).LT.TOL)GO TO 14 

X=X1 

GO TO 15 
14 P 4P 4 I D = F 2 
RHW CR 4 = F 3 
W 4WC 4=X 


W CR 4 = WCK 3* SQR T ( T4T3) 75 

W4 = WCR 4*W4WC4 

VL'4 = U4-W4*ABS( S I N I F 1* B 4 ) ) 76 

WU4=W4*SINI FI* B 4 ) 77 

T 4T 4=1. + G 3* ( U4**2+2.*U4*W U4 ) / WCR4* *2 

P4PPDP=P 1PP0P*P4P4ID*I TT3*T4T3*T4T4>** (G/G2)*P2PP1 P 78 

SINB4=SIN<Fl*B4) 79 

CrSB4=CHS( F1*B4) 80 

VCR4 = WCR4*SQRT( T4T4) 81 

V4 = S(0RT ( W 4**2* CO SB 4** 2 + ( U4+W4*SI N B4 ) * *2 ) 82 

P 4P 4P = ( 1 .-G3*{ V4/VCR4)**2 )**(G/G2 i 83 

P4P0P=P4PP0P*P4P4P 


6C 
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T4P=T4T4*T4T3*T3PP 
TEFTCR=( VCR0/IC19.46)**2 
WK— VIJ3*U 3— VU4* U4 
W K = to K / AC/XJ 
WK1=CP*< TO— T 4P ) 

W I DP =CP * TO* ( l.-P4PP0P**(G2/G) ) 84 

W ID=CP*TO*l 1 P 4PUP** (G2/G )) 85 

ETAT=toK/WlDP 
ETAS = /JK/WID 

P 3P4P = ( 1 G 3* ( V 3VC 3**2— (VIJVC3**2— { U3/ VCKO) **2 ) **2) )**( G/G2 ) 86 

P3P3P=( l.-G3*V3VC3**2)**<G/G2) 87 

P4PP3P=P4PPOP/P2PPOP 

WKI0P = CP*T0*{ 1.-P4PP3P** (G2/G) )- ( U3-VU3 ) **2/ < 2 . * AC*X J ) 88 


RrTJ.OS=XM*XK*t W3**2+W4**2> / (2.*AC*XJ) 

vc-vc*o*vnvcn 

V 2= VCRO* V2VC2 

STALOS=XK*( VO**2+V2**2) /2. /AC/XJ 
C ER =i<0 TLO S /W I D 
CES=STALOS/toID 
DEL=PO/2116.22 

W EO = Sail* SORT ( THE TCR ) /DEL*E P S 89 

W JEQ=frfK/THFTCR 

DEE = V4**2/( 2.*AC*XJ*ViID) 

E1AF=ETAS*0FF 
XI F=V4** 2/2 ./XJ/AC 

OEI=(J3*SIN(8 3-PHI ) )** 2/ ( 2. *AC*X J*WI 0 ) 90 

XL I=DEI*WID 

A L 2= 5 f . 3*AL 2 

R3=57. 3*83 

P 1-1 = 57 .3*PHI 

FTAl = FTAE +0E I 

E T AR =£ TA I +DEE 

XN=60.*U3/PI/D3 


FCSPD=XN/SQRT1 THETCR ) 91 

XNU = IJ3/S0RT <?.*AC*Xj*toID> 92 

SPP AR = SW*XN*FP S/DE L 
P FP P4P= 1 ./P4PP0P 

X NnT = XN / SOR T ( TO) 93 

W T O P = S W * SO R T ( TO ) * 1 44 • /PO 94 


XT0P=SW*WK*3. 5616* 1 44 . /XN /PO*AC*XJ 
P 4= P4P0P *P0 

T 4 = T4P* ( 1.-G3* (V4/VCR4)**2) 

RF0A=P4/R/T4 

04=SW/RH04 

Cr = lJ3/S0RTC 2.* AC*XJ*toIDP> 95 

H A = X J* W I DP 

XN SP = XN* SOR T(Q4)/H4**.75 96 97 

SPP AR=Sw*XN/DEL 

G00 = 7537 • 129*WK*Sto/PO/XN*AC*XJ 
POP P 4=1. /P4P0P 
I F< MODE .NE .2 ) GO TO 20 
EX=FTAS 

IFt ETESr.FQ.2) EX=E TAT 
IFIK.E0.6) GO TO 25 

CALI GFTMi XK»XMf2tEX*K) 106 

GO TO 23 
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u u 


2 1 V lVCl = ViVCl40ELV 

I FI VlVCi .GT .VEND ) GO TO 22 
GO TO. 23 
22 V1VCI=VSTART 
DFL V=OEL V l 

I F ( MOi)E .60*1) GO TO 5 
I F( MODE «EQ *2 ) GO TO 3 
I F< MOOE.EO .3 i GO TO 4 
2 5 CONTINUE 


WRITE! 6, 

5 C 9 ) 

XM , 

XK , XK 1 , PE RC T , 

U 3 ,U4 , PREF , WREF 

123 

IF! PTEST 

• EO, 

1 ) 

WR I TE ! 6,529) 

PDSGN 

124 

IFtPTFST 

•EG). 

2) 

WR ITE ( 6,526) 

PDSGN 

126 

I F! ETFST 

• EO. 

1) 

WRITE (6,527) 

ETAD 

128 

IF! FTFST 

• EO. 

2 ) 

WRITE (6, 528) 

ETAD 

130 


GO TO 20 
C CHDKFC 

13 I F ( MODE • EO • 3 ) GO TO 20 
V IV 0 1 = V1VC 1-DE L V 

0 EL V=OEL V/ 2 • 

1 F ( DFLV.LT.5.F-4) GO TO (5,3), MODE 
GO TO 21 


OUTPUT FOR SPEED V1VC1 

20 CON T l MUF 

WRITF(6,505) 

505 FORMAT { 1 F9 ) 

WRI TF( 6, 506 ) V 1 VC i , P 2P P OP »V3,WCR3,VCR0, P2P PI Pt 03 ,W4 ,RHOOP, 

1AL 2, TT3, WCR4,PVCR0,R2, T3PP , VU4 ,RHD VC i , VUVC 3 , T4T3 ,T4T4, SW, 

?VXVC3,W 3 WC 3 » P4PP OP , RHO VCQ • V3VC3*RHXVC3t VCR4 

WR TTE(6, 50 7 ) VO VCO , VU3 , P4P4I D, V4. V2 VC2 , W3 , RHWCR4 , P4 P4 P, 

1RF2VC2, VX3,W4WC4,P4P0P,WK,ETAE , WI 0 P , T4P , WI 0 • ET AT ,FTAS , WEQ, 

2ETA T, XL F *OEF , XL I #ROT LOS , S TALOS , DE R , DE S , XN , EOSP D ,XNU, THFTC R 

W R I T E ( 6 • 50 8 ) DEL , SPPAR , W1E C ,P0PP4P , XNOT , WTOP ,XTOP,DEI . 

IP OPP 4 , PH I » P 1PP0P.LTAR ,XNSP 

GO TO 21 
C OUTPUT INPUT 

40 WRITE! 6. 501 ) AO , A1 • A3 ,A4,D1»D3 , XL ,AL1*B4,ZBR * PO,TO, 

1 U 3 1 , IJ4 1 • G * R ,CP * E P S 

r F! MODE •FO . 1 ) WRITE < 6, 5C2) 

I F { MOt)E . FQ * 2 ) WRITE 16,503) 

I F ( MOD t * EO • 3 ) GO TO 525 

509 FORMAT! 1 4X , 1 HM * G 15 .6/ 1 4X , 1HK *G 1 5 .6/13X , 2HK1 , G1 5. 6/ 7X ,8H PERCENT, 
1G15.6 / nX,2HU3,GL5.6/13X,2HU4,G15.6/LX,14HDESIGN PI ,/PO, G15.6 / 

2 IX , 14F0ES IGN SMALL- W. Gl 5, 6 ) 

WR ITF! 6, 509 ) XM,XK,XK1 , PE RC T , U 3 , U4 , PREF , WREF 
WRITE! 6,505) 

Gf TO 50 

525 WR I TE! 6, 504 ) 

GC TO 50 

C FORMATS 

529 FORMAT ( 2X, ) 3HDESIGN P0,/P4 G15.6) 

526 FORMAT! IX, 14HDESIGN P0,/P4, G15.6) 

527 FORMAT! 2X 1 3 FDE SIGN ETA-S G15.6) 


140 


141 


142 


143 


14 5 
146 
148 


152 

153 

155 
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528 FORMAT { 2X 13HOESIGN ETA-T G15.6) 

400 FORMATS BE 10.5) 

501 FORMAT < IF 1 . 

l 36 X * 49HQ FF DESIGN PERFORMANCE OF A RADIAL INFLOW TURBINE // 

12JH TURBINE GEOMETRY AO , G 16. 6 ,2X , 2 HA 1 , G16 . 6, 2 X, 2H A3 ,G16 . 6 , 2X , 
22FA4,G16.6. 2X. 2HD1 ,G16.6 / I9X .2HD3.G16. 6 ,2X ,1HL ,G 17 .6 , 2X , 3 HAL l • 

3G 16*6* IX .2HB4.G16.6.2X .IH 2,016.6 /21H INLET CONDITIONS PO, 

4G 16 • 6. 2X.2H TO, 016*6 /14H DESIGN SPEEDS ,5 X , 2HU3 ,G16 .6 ,2X , 2HU4, 
5G16.6/20F FLUID PROPERTIES G , G 1 7. 6 , 2 X ,1 HR * G 17 .6 , 2 X f 2HCP ♦ G 1 6 . 6 , 

6 2 X * 3FFPS.G16.6) 

401 FORMAT { I 1 ) 

502 FORMAT I 8HKM0DE 1* 2X , 3 7HPERF0RMANC E RUN VARYING PERCENT SPEED ) 

503 FORM AT ( 6HKM0DE 2, 2X , 3 IH MAN UAL DETERMINATION OF M AND K ) 

504 FORMAT ( 8HK MODE 3, 2X , 34HALTCMA TIC DETERMINATION OF M AND K ) 

506 FORM A T( 12H V/VCR ) 1 ,G 1 5. 6, 9 X , 7H P2 , /PO , , Gl 5 . 6 , 14 X , 2HV3 *G 1 5 . 6 , 

l l lX*4HW0R3tG15.6/7X*5HVCR)0»G15.6*9X*7HP2 ,/Pl . ,G15 .6 , 1 IX, 5HBET A3 , 
2G15.6. 13X #2 HW 4tGl5.fi/8X* 5H RHCU CNGi5.6*lOX* 6HAL PHA2 ,Gl5.6*8X f 
38 FT* , / T* ) 3.G15. 6. 11X.4HWCR4 ,G15.6/ 12H RHO-VCR) 0, G15.6 , 14X , 2HR2, 
4G15.6. I 2 X , 4HT 3 » * .Gl 5.6 * 12 X ,3HV U4 ,GL5 .6/12H PV/ PVCR ) 1 , G15 .6* 

58 X, BHVU/VCR ) 3 , G l 5. 6 . 7X . 9H T 4 , ,/T3* . ,G l 5 . 6 , 7 X , 8HT4 ,/T4 f , ,G15.6 / 

61 2H SMALL-W.G15.6.8X, 8HVX/ VCR) 3 ,G15. 6 , 9X, 7HW/WCR 13 ,G15 .6 , 6X , 

77FP4./PO. .G15.6 /12H P V/P VCR ) 0 .G 1 5 . 6 , 9 X * 7H V/ VC R ) 3 , Gl 5 . 6 , 6X * 

8 1CHPVX/PVCR )3*G15.6* ILX.4HVCR4 ,G1 5 • 6 ) 

507 FORM A T ( 1 2H V/VCR ) 0.G15 .6 , 1 3 X , 3HVU3 V G 15.6 ,9X , 7HP4/P4ID, G15. 6* 

1 1 t 2HV4 * G 15.6 /12H V/VCR) 2 *Gi 5 . 6 , 14 X , 2HW3 ,G15.6,7X, 

29FPW/PWCR )4,G15.6,9X,6HP4/P4,,G15. 6 /12H PV/ PVCR ) 2 ,G 15 .6 , 13X , 

3 3 F V X 3 , G15.6.9X, 7HW/WCR ) 4 *G1 5. 6 ,9X,6HP4/PC, .G15 .6/ 1 OX ,2HW 1 • G15 .6* 

41 IX. 5HET A—E.G15.6, 12 X . 4HW . I D .G 15. 6 . 1 2X , 3HT4 . .G15 .6/9X, 3HWID.G15.6, 
51 IX. 5HETA-T.G15.6, 1 IX.5HE TA-S, G15. 6 *1 2 X , 3HWE0, G15 . 6 / 

61 F 6X. 5HFTA-I.G15. 6 , 14X , 2HLE , G 15. 6 ,13 X ,3 HDEE , G15 .6 , 13X . 

72 FL 1,0,15.6/ 12H ROTOR-LOSS ,G1 5 . 6,5 X . 1 IH S TATOR- LOSS ,G15 .6, 

81 2X, 3FDER, G15.6, 12X.3HDES.G15. 6/1H ,6X, 5HSPEED ,G15 .6 ,8X, 

X8FE0— SP FED, 

9G15.6. 14X.2FNU.G15.6.7X.8H THE T A— CR ,G15.6 ) 

508 FORMAT! 7X , 5HDELTA , G 15 . 6 , 7X . 9H SPEED-PAR .G15.6 , 

11 IX * 5HW 1— EQ *G 1 5.6* 8X.7HP0 ,/P4, ,G15 .6/9X, 3HN0T, G15.6 ,12X , 4HWT0P, 
2C15.6. 13X.3FTOP.GI5.6.12X. 3HDE I ,015.6/ 

31 F 5X • 6HP0, /P4*Gi5.6*13X* 3 HP HI ,G 15 . 6 , 9 X , 7H P 1 , / PO , , Gl 5 . 6 , 

4 1 OX, 5HETA-R, G15.6/ IH 8X , 3HNSP ,G15. 6) 

END 


Internal Program Variables and Definitions 


AO 

A 0 

A1 

A 1 

A3 

A 3 

A4 

A 4 

AC 

acceleration constant 

AL1 

a x 

AL2 

“2 

B3 

^3 

B4 

^4 

COSAL1 

cos 

COSB3 

COS /3g 

COSB4 

cos £ 4 

CP 

Cp (input variable) 

D1 

D 1 

D3 

D 3 

DEE 


DEI 


DEL 

6 

DELV 

incremental value of (v/v\ 1 (input variable) 



DER 

A % 

DES 

<1 

DM 

Am 

ETAD 

design value of efficiency 

ETAE 

% 

ETAI 


ETAR 

% 

ETAS 

% 

ETAT 


EPS 

e 

F 

temporary storage 

FI 

temporary storage 

F2 

temporary storage 

F2P 

temporary storage 

F3 

temporary storage 

F3P 

temporary storage 

FP 

temporary storage 

G 

r 

G1 

r + l 

G2 

r - i 

G3 

(r - 1)/ (r + !) 

H4 

temporary storage 
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KMAX 

MODE 

PO 

POPP4 

POPP4P 

P1PPOP 

P2PPOP 

P2PP1P 

P3P3P 

P3P4P 

P4 

P4POP 

P4P4ID 

P4P4P 

P4PPOP 

P4PP3P 

PDSGN 

PERCT 

PHI 


maximum number of k's 

if MODE = 1, performance run varying percent speed 

if MODE = 2, manual determination of m and k 

if MODE = 3, automatic determination of m and k (input variables) 

PO 

P(/?4 

P(/P4 

p’i/po 

p^p’o 

Pg/Pi 

P 3 /P 3 

P 3 /P 4 

P 4 

P 4 / p 0 

P4/P’4,id 

P 4 /P 4 

P 4 /P 0 

P 4 /P 3 

specified pjy/p'^ or p'q/p 4 (input variable) 
percent of design speed (input variable) 

<P 
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PIN 

PREF 

PVCRO 

Q 

Q1 

Q2 

Q3 

Q4 

R 

R1 

R2 

R3 

RH2VC2 

RHOOP 

RHOVCO 

RHWCR4 

RHXVC3 

ROT LOS 

SINAL1 

SINAL2 

SINB3 


pJ/144 

estimate of Pj/Pq to determine (input variable) 

(P’ v cr)o 

temporary storage 
temporary storage 
temporary storage 
temporary storage 
temporary storage 
R (input variable) 

r l 

r 2 

r 3 

(pV/p’V cr ) 2 

p 0 

(pV/p'Vc^o 

(pw/p"w cr ) 4 

(pv x /p'v cr ) 3 

l e 

sin 0*2 
sin a 2 
sin /3g 
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SPPAR 

wN /6 

STALOS 

L S 

TO 

T ’o 

T4 

T 4 

T4P 

T 4 

T3PP 

rpf T 

1 3 

T4T3 

T 4/T 3 

T4T4 

(t 7 t ") 4 

THETCR 

e cr 

TT3 

(t' 7 t ’) 3 

U3 

U 3 (input variable) 

U31 

initial U 3 

U4 

U 4 (input variable) 

U41 

initial U 4 

VO 

V 0 

V2 

V 2 

V3 

V 3 

V4 

V 4 

vovco 

( V / V cr)o 

V1VC1 

see input 



V2VC2 

( V / V cr) 2 

V3VC3 

( v / v cr)s 

VCRO 

V cr,0 

VCR4 

V cr,4 

VEND 

end value of ^V/V 

VSTART 

starting value of | 

VU3 

V u,3 

VU4 

V u,4 

VU30PT 

eq. (B16b) 

VTJVC3 

(V V cr ) 3 

VXVC3 

(V V cr) S 

W3 

W 3 

W4 

W 4 

W3WC3 

( W / W cr) 3 

WCR3 

W cr,3 

WEQ 

W eq 

WEQD 

see input 

WID 

Ah ld,s 

WIDP 

Ah ld,t 

WREF 

estimate of W tc 
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WTOP 

wyn/p’o 

WU30PT 

eq. (B16c) 

WU4 

W u, 4 

WX3 

W x, 3 

X 

temporary storage 

XI 

temporary storage 

XJ 

constant, 778. 16 

XK1 

k l 

XLE 

l e 

XN 

N 

XNOT 

n/t'q 

XNSP 

NSP 

XNV 

V 

XTOP 

r/p! 

' F in 

ZBR 

Z (input variable) 


Subroutine GETK1(XK1:X1) 

This routine evaluates the loss coefficient and a design (v/ V. cr ^ ^ a 

specified Pq/Pi and weight flow w where 

XK1 loss coefficient kj 
XI design (v/V cr ) l 
The subroutine is as follows: 
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S IBR01JT INF GFTKHXK1.XIJ 

COMMON / GETK/ P , W ,G . AO , A 1 . AL , T . R, PO 

SUB( X, G ) *X*( l.- ( G- 1 . )/ (G+l . )*X*X) ** ( 1 . /( G-l. ) ) 

VCRO = ( 64.34*G*R*T /(G+l.))**. 5 

RECO=PO/R/T 

Q=W/(RHOO*AG*VCRO) 

ROCO=Q 
X = . 1 

1 F = SUB( X, G J-Q 
F 1 =( F+O ) /X 

FP=-?.*X*X/(G+1. )*F1**(2.-G)+F 1 
X 1 = X— F/FP 

I F ( ARSI (XI— X)/X1).LT. l.E— 5 ) GO TO 2 

X=X1 

GG TO 1 

2 vr=xi 

Q=W/(P*KHOO*VCRO*A1*COS( AL*. 01 7451 ) 

X = .l 

3 F = S(JB ( X, G )— 0 
RDC1=F+Q 

F 1= { F+ 0 ) / X 

FP=-2.*X*X/< G+l . )*F 1**( 2.-G1+F 1 
X 1=X— F/FP 

IFIABSH X1-X1/X1 l.LT. l.E-5) GO TO 4 
X = X1 
GC TO 3 

4 XK1 = (P**(CG-1.)/G)-1.1*(1.-(G-1.)/(G+1.) *X1**2 >/ 
1( ( G— 1 . ) / ( G + 1 • ) * ( VO** 2+Xl**2 ) ) * (-1. ) 

RETURN 

FND 


Subroutine GETIVMXK, XM,IGO, ETEST, K) 


For IGO = 1, this routine computes a value of XM as a function of XK by using 
equations (B17b) to (B24b). For IGO = 2 this routine varies the value of XK by false 
positioning until all design specifications are met. The following are defined for this 


subroutine: 


XK 

loss coefficient 

XM 

loss coefficient 

IGO 

index for mode of operation 

ETEST 

design value of or jj s 

K 

indicator used in method of false positioning 


The subroutine is as follows: 


32 


S18R0UT INE GETM (XK.XM , IGO.E TEST,K) 

COMMON /GETK /PR EF ,WREF ,G »A 0 , A 1 , AL1 , TO, K, PO, P4PPOZ ,CP ,G1 , G2, G3, 
1ETAT,U3,IJ4,VU3,B4, T4 T 3 , fcC R3 ,TT3 ,W3V*C3 , B3 *PHI , P2 P PI P , PI PPOP 


I F ( IGO.GT. 1 ) GO TO 100 
P4PPDP= 1 ./P4PP0Z 

DPI D=CP*TU* ( 1 .-P4PP0P**(G2/G) ) 5 

DE = CMO*ETAT 

V l,4 = < OK 2*50 3 3 . - V U3*U3 ) /U4 

V L4=-VU4 

W4=(VU4— U4) /AB S ( SI N ( • 0 1 74 5*B4 ) ) 6 

U 4=-W4 

V* CR4 = irtlCR 3*SQRT{ T4T3 ) 7 

W4WC4=W4/uJCR4 

fcL4=H4*S IN ( .01745*8 4) 8 

T 4 1 4= 1 . +G 3*{U4**2+2.*U4*WU4 ) /WCR4**2 

P4P4I F) = P 4PP0P/ ( P2PP 1P*P 1 PPOP* ( TT3*T4T3*T4T4)** (G/G2) ) 9 

T=W3WC3**2/T4T3 

T ]=-( P4P4 IC**( G2/G )-!.)*( 1 .-G3 *W4WC4**2 ) /G3 10 

XMK=< Tl-T*( SIN (B3-PHI ) )**2) / ( W4WC4**2+ T ) 11 


XM=XMK/XK 

RETURN 

100 GG TO < IC1. 102) .K 

101 X K 1 — XK 

E 1 = F TE ST 
D IF1=ETFST-FTAT 
K =2 

XK=XK+.005 

RETURN 

102 X K 2= XK 

F 2=ETF ST 

C IF2 =E TEST- ETA T 

I F ( D I F2* D IF 1 ) 106,105,104 

104 E 1 = E 2 

X K 1=XK 2 
X K = XK + .00 5 
D IF 1=0 IF 2 
RETURN 

105 XK=XK2 
K =6 
RETURN 

106 XK=- ( XK2-XK 1 )/ { D IF2-DI F 1 ) *DIF l + XKl 
K =6 

RETURN 

FND 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, December 6, 1968, 
125-23-02-10-22. 
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APPENDIX A 


SYMBOLS 


A 

C I 

D 


g 

Ah’ 


Ah ld,t 


J 

kp k, m 


L 

l 

N 


NSP 



r 

T 

U 

V 

W 


w 

Z 


a 


area, sq ft 

specific heat at constant pressure, Btu/(Lb)(°R) 
diameter, ft 

o 

dimensional constant, 32. 1741 ft/sec 
specific turbine work, Btu/lb 

ideal turbine work based on inlet total to exit static pressure ratio, Btu/lb 

ideal turbine work based on inlet total to exit total pressure ratio, Btu/lb 

mechanical equivalent of heat, 778. 029 ft-lb/Btu 

loss coefficients, dimensionless 

kinetic energy loss, Btu/lb 

axial distance, ft 

turbine speed, rpm 

specific speed 

pressure, lb/ sq ft 

turbine inlet pressure, psia 

gas constant, (ft -lb)/ (lb)(°R) 

radius, ft 

absolute temperature, °R 

blade speed, ft/sec 

absolute velocity of gas, ft/ sec 

gas velocity relative to rotor, ft/sec 

mass flow rate, lb/ sec 

number of rotor blades at tip of rotor 

absolute gas angle, angle between absolute velocity vector and meridional 
plane at mean channel, positive when direction of tangential velocity 
component and direction of wheel velocity agree and negative when they 
disagree 
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I 


0 

r 

r 

6 

€ 


Ve’Iv’Ir 


% 

\ 

6 ny% 

cr 

v 

<P 

P 

Subscripts: 

cr 

E 

eq 

I 

id 

opt 

R 

S 


relative gas angle, angle between the velocity vector relative to wheel 
and meridional plane at mean channel; same sign convention applies 
as for a 

torque, ft -lb 

ratio of specific heat at constant pressure to specific heat at constant 
volume 

ratio of turbine inlet total pressure to U.S. standard atmospheric 
pressure, pJj/2116.22 

function of y used in relating parameters to those using air inlet 
conditions at U.S. standard sea-level conditions 

ll y/ (r ' l) j( y* + l^ j y * /(r * _1) ~ 

modified efficiency figures used in studying loss distribution in turbine, 
dimensionless 


efficiency based on ratio of inlet total to exit static pressure 

efficiency based on ratio of inlet total to exit total pressure 

squared ratio of critical velocity at turbine inlet to critical velocity at 
U.S. standard atmospheric temperature, ^V cr ^ 1019. 4(^ 

blade -jet speed ratio, Ug^/2gJ Ah!^ g 

optimum rotor incidence angle, deg 

gas density, lb/ cu ft 


condition corresponding to Mach 1 
exit 

air equivalent (U.S. standard sea level) 

incidence 

ideal 

optimum 

rotor 

stator 


u tangential component, positive when its direction agrees with that of 

wheel velocity and negative when it disagrees 
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X 


meridional component, component in plane containing axis of rotation 

0 station at turbine inlet 

1 station immediately upstream of stator exit 

2 station where flow from stator exit is assumed to occupy entire cylindrical area 

3 station immediately upstream of rotor inlet 

4 station immediately downstream of rotor exit 
Superscripts: 

(') absolute total state 

(") total state relative to rotor 

(*) U.S. standard sea-level conditions (temperature, 518. 67° R; pressure, 

14. 696 psia) 
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APPENDIX B 


THERMODYNAMIC EQUATIONS 

The equations in this appendix are listed in the order used by the computer. A 
complete set of calculations is made for each specified value of (v/ V cr ) j at each speed. 


pV \ = / V \ 

)’V / J 

ci/ ^ ' cr/ ^ 


y - 1 /v 

y + 1 IV 


cr/ 


l/(y-i) 


(Bl) 


Equations (B2) to (B5) are to be solved simultaneously for Pj/Pq: 


V 


p 0 \p’ v , 


cr/ 


( P ’ V 4 


cos 


(B2) 


PV 


w 


^p’Vcr/o (P’ V cr A >0 


(B3) 


Calculate 


( v / y cr)o 


from 


r PV 
P’V, 


cr/ 


cr/ 


y - 1/ V 
y + i \v ( 


cr/ 


1/ (y-l) 


(B4) 






p 0 


y-l 

y + 1 


v\ 2 + / v^ 2 


' V cr/ 0 V^cr/ 1 


y/ (y-l) 




y - if v 

y + 1 \V 


cr/ 






(B5) 
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(B6a) 



P2 

P 2 = Jl®. 

p’l Pi 

Pi) 



sin Q!j 

Yo = Y 1 

sin tig 


(B6b) 


(B7) 


(B8) 


(B9) 


(BIO) 
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(BH) 



(B12) 


(B13a) 


(B13b) 


(B14a) 


(B14b) 


(B14c) 


(B14d) 
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Calculate the angle 



as follows: 

V u, 3, opt _ 1 _ 1. 98 

u 3 z 


(B14e) 


(B14f) 


(B15a) 


0315b) 


(B15c) 


(B15d) 


(B15e) 


(B16a) 
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V u, 3, opt “ U 3 


(B16b) 


V u, 3, opt\ 

, U 3 / 

W u, 3, opt = V u, 3, opt “ U 3 

y = tan- lWu 7^°P t - 
V x,3 

The following three equations are to be solved simultaneously for (W/W ) 4 and 

r t / it \ / ’ 

P/P 4 > l d : 




(B16c) 

(B16d) 

(B17a) 


(B17b) 

(B17c) 
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(B18a) 



W u 4 = sin /3 4 (j3^ will be negative) 


P 4 PgPi P 4 


PO PlP()P4,id 



y/ (y-i) 


1 cx, 4 ^cr, 4 



A 1/2 


w 


4 


= W 


cr, 4 



4 



P 4 

P4 


y - l 

y + 1 



y/ (y-l) 


(B18b) 

(B18c) 


(B18d) 
(B18e) 
(B 19a) 


(B19b) 

(B19c) 


(B19d) 


(B20a) 



(B20b) 


P 4 = P 4 P 4 

PO PqP 4 


V u.4 = U 4 + W 4 sin ^4 


rpf rpT t 

T r _ 1 4 1 4 T *» 
A 4 rp" ip" 3 
1 4 1 3 


... \.3<V V u,4 U 4 

gJ 


Ah id, t “ c p T 0 


1 - 


v p oy 


(y-l)/y 


Ah id, S _ c p T 0 


i-i* 


(y-l)/y 




Ah’ 


Ah. 


Va = 


id, t 
Ah' 


s Ah! 


id, s 


(B21a) 

(B21b) 

(B22) 

(B23a) 

(B23b) 

(B24a) 

(B24b) 


The foregoing equations have given the flow conditions through the turbine and the 
efficiencies. The following equations give quantities which are useful for special pur- 
poses, such as identifying sources of loss in the turbine: 


L r = mk 


W^ + V/\ 

2gJ 


(rotor loss) 


(B25) 
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V 


(B26a) 


V 0 V cr, 0 . v 


CT 


V 2= V cr,0. v 


cry 


V 0 + V 2 

L q = k — - (stator loss) 

S 2gJ 


A % = 


Ah! 


B — (fractional rotor loss) 

id, s 


Ajjg = 


Ah! 


(fractional stator loss) 


id, s 


L™ = — - (exit loss) 

E 2gJ 




Ah! 


_I! — (fractional exit loss) 

id, s 


Wo sin^ (/3„ - <p) 

L t = - (incidence loss) 

1 2gJ 




Ah 


(fractional incidence loss) 


id, s 


(B26b) 

(B26c) 

(B27) 

(B28) 

(B29a) 

(B29b) 

(B30a) 

(B30b) 
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n E = v s + ±v E 


~\ 


Vr- Ve + A % > 


r tl = VR + A Vl 


e =\ cr >° 
cr \1019. 46 j 


e = 


1.40 


W) 


y/ (r-l) 


1.8929 


6 = 


'0 


2116.22 


w eq = w V. c £. e (equivalent weight flow) 


—— e (weight flow-speed parameter) 
6 


N = — — — (equivalent speed) 

6,1 /v 


U, 


v = 




(blade -jet speed ratio) 


Ah’ = — - (equivalent specific work) 
e q q 

cr 


NOT = 


N 


7 % 


(B31) 

(B32a) 

(B32b) 

(B32c) 

(B33) 

(B34) 

(B35) 

(B36) 

(B37) 

(B38) 
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wJt! 


WTOP = 


0 


Pi 


in 


(B39) 


where p! n = Pg/144 


TOP = 


where r = 89 155 w Ah’/N (in. -lb) 


m 


(B40) 


NSP = 


H 


3/4 


(specific speed) 


(B41) 


where Q = w/p^ and H = J Ah.^^.. 


For calculation of 


P 4 = 



T 


4 



1 


y - 1 

y + l 



P 4 = 


P 4 


RT 


4 
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